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1 INTRODUCTION 

In accelerators with very short bunches, such as is envi- 
sioned in the undulator region of the Linac Coherent Light 
Source (LCLS)[[l]], the wakefield due to the roughness of 
the beam-tube walls can have important implications on the 
required smoothness and minimum radius allowed for the 
beam tube. Of two theories of roughness impedance, one 
yields an almost purely inductive impedance^, the other a 
single resonator impedance [^J|; for smooth bunches, whose 
length is large compared to the wall perturbation size, these 
two models give comparable results^. 

Using very detailed, time-domain simulations it was 
found in Ref. [Q] that a beam tube with a random, rough 
surface has an impedance that is similar to that of one with 
small, periodic corrugations. It was further found that the 
wake was similar to that of a thin dielectric layer (with 
dielectric constant e sw 2) on a metallic tube: W z (s) w 
2/Cq cos kos, with wave number and loss factor 



fcn = 



and /Cn 



Z c 
2na 2 



(1) 



with a the tube radius, S depth of corrugation, and Zq = 
377 ft. For the periodic corrugation problem this result 
was inferred from simulations for which the period p ~ <5. 
On the other hand, at the extreme of a tube with shallow 
oscillations, with p ^> <5, the impedance was found, by a 
perturbation calculation of Papiernik, to be composed of 
many weak, closely spaced modes beginning just above pi 
phase advance[|[. 

In this report we find the impedance for two geometries 
of periodic, shallow corrugations: one, with rectangular 
corrugations using a field matching approach, the other, 
with smoothly varying oscillations using a more classical 
perturbation approach. In addition, we explore how these 
results change character as the period-to-depth of the wall 
undulation increases, and then compare the results of the 
two methods. 

2 RECTANGULAR CORRUGATIONS 

Let us consider a cylindrically-symmetric beam tube with 
the geometry shown in Fig. [j]. We limit consideration here 
to the case 5/a small; for the moment, in addition, let 
5/p > 1. We follow the formalism of the field matching 
program TRANSVRS[f|: In the two regions, r < a (the 
tube region, Region I) and r > a (the cavity region, Region 
II) the Hertz vectors are expanded in a complete, orthogo- 
nal set; E z and are matched at r = a; using orthogonal- 
ity properties an infinite dimensional, homogeneous matrix 
equation is generated; this matrix is truncated; and finally, 
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Figure 1 : The geometry considered. 

the eigenfrequencies are found by setting its determinant 
to zero. We demonstrate below that, for our parameter 
regime, the system matrix can be reduced to dimension 1, 
and the results become quite simple. 

In the tube region, the z -component of the Hertz vector 



with Iq the modified Bessel function of the first kind, and 

2nn 
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p 



X, 



(3) 



with fio the phase advance and k the wave number of the 
mode. In the cavity region, 
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cos[a s (z + g/2)}, (4) 
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with Kq the modified Bessel Function of the second kind. 
E z and are given by 



E, 



dz 
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-jk-^— ■ (7) 
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Matching these fields at r = a, and using the orthogo- 
nality of er^ nZ on [— p/2,p/2], and cos[a s (,z + g/2)] on 
[— g/2, g/2] we obtain a homogeneous matrix equation. To 
find the frequencies, the determinant is set to zero; i.e. 



det 
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with the matrix N given by 



sin(/3„«?/2) 

{PI ~ arftg \ cos(/3„<?/2) 

and the diagonal matrices 1Z and X by 

R'r 
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s odd 
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(9) 



TZ S = (l+5 sa )ka 



X n = kc 
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For the beam, on average, to interact with a mode, one 
space harmonic of the mode must be synchronous. We will 
pick the n = space harmonic to be the synchronous one; 
i.e. let /?o = k (we take the particle velocity to be u = 
c). We truncate the system matrix to dimension 1, keeping 
only the n = and s = terms in the calculation (the 
other terms are small). Now if kS is small, then the s = 
term in 1Z becomes TZg — 2/(k5), the n — term in I is 
2q = ka/2, and Nqq ~ 1. Eq. 0then yields 



k = 



2p_ 

aSg 
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which, for p = 2g, equals fco of Eq. [j]. 

The loss factor is given by K = \V\ 2 /[4Up(l - /? 9 )][0, 
with V the voltage lost by the beam to the mode, U the 
energy stored in the mode, and j3 g the group velocity 
over c. The voltage lost in one cell is given by the syn- 
chronous (n = 0) space harmonic: V = A^p, and the 
energy stored in one cell, U = 1/(2Zqc) j E ■ E* dv, 
is approximately that which is in the n = space har- 
monic: U = nAla 2 p{l + k 2 a 2 /8)J(2Z c) (for details, 
see Ref. [||]). For (3 g , we take Eq. || truncated to dimen- 
sion 1, and expand near the synchronous point. Taking the 
derivative with respect to /3q and then setting [3q — k we 
obtain: 



(12) 



The loss factor becomes K- = K-q. 

The above method can be extended to modes of higher 
multipole moment m, in which case the beam will excite 
hybrid modes rather than the pure TM modes of above [jj. 
Again the system matrix can be reduced to the n = and 
s = terms, and the lowest mode wave number and loss 
factor have a simple form (for 1 < to -C a/S): 



k= /fm + l)p and K 
V ad 9 



Z Q c 
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(13) 



and (1— j3 g ) — m(m + 2)5g/ (ap). We note that the dipole 
(to = 1) frequency is equal to the monopole (m = 0) fre- 
quency. Also, the wake at the origin is the same as for the 
resistive-wall wake of a cylindrical tube[||], as we expect. 

Running TRANSVRS with a matrix of dimension 40, 
we obtain a typical dispersion curve (see Fig. |2|). Here 
k/ko = 1.07, K./K.Q = .94. Note that even when S/a is 
not so small, e.g. for bellows with 8/a ~ .2[0], the analyti- 
cal formulas are still useful. Fig. |5] shows how the strength 
and frequency of the mode change as the period of undula- 
tion is increased. The scale over which JC drops to zero is 
Po w 7r yj a5g j 2p. By p ~ po, the one dominant mode has 
disappeared, and we are left with the many weak, closely 
spaced modes, beginning just above kp = ir, that were 
found by Papiernik. 

3 SINUSOIDAL CORRUGATIONS 

Let us assume now that the pipe surface is given by 
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Figure 2: Dispersion curve example. 
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Figure 3: An example showing the effect of varying p. 



where 2tt/k is the period of corrugation, and h is its am- 
plitude. We assume that both the amplitude and the wave- 
length are small, h <C a and na 3> 1. This allows us 
to neglect the curvature effects and to consider the surface 
locally as a plane one. We will also assume a shallow cor- 
rugation Iik <C 1, i.e. the amplitude of oscillation is much 
smaller than the period. 

Introducing a local Cartesian coordinate system x, y, z 
with y = a — r (directed from the wall toward the beam 
axis), and x directed along <j>, the surface equation becomes 
y = y (z) = h sin kz. The magnetic field near the surface 
H x (y, z) does not depend on x (that is <f>) due to the ax- 
isymmetry of the problem. It satisfies the Helmholtz equa- 
tion 

*H. .*H..,.„ (15) 



dy 2 dz 2 
with the boundary condition 



k 2 H r 







(nVH)\ y=ya = 



(16) 



where n is the normal vector to the surface, ft = 

(0, 1, — JlK COS Kz). 

Note that the longitudinal electric field E z can be ex- 
pressed in terms of H x , 



= a — h sin / 



(14) 



i dH x 
E z = --- 
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(17) 



Using the small parameter h/a, we will develop a perturba- 
tion theory for calculation of H x near the surface and find 
how E z is related to H x . 

In the zeroth approximation, the z dependence of H x is 
dictated by the beam current periodicity, 



H x (y,z)=H(y)e 



i k z 



(18) 



Putting Eq. ^ into Eq. ([B|) we find that d 2 H/dy 2 = 0, 
hence Ti(y) = Ho + Ay, where the constant A can be 
related, through Eq. (|T7|), to the electric field on the surface, 
A = ikE z . We will see below that A is second order in h. 

For a flat surface, for which n = (0,1,0), from the 
boundary condition (16), we would conclude that A = 0, 
however, the corrugations result in a nonzero A, and hence 
E z . Substituting the magnetic field (|l8|) into the right hand 
side of Eq. © one finds 
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Clearly, the boundary condition is not satisfied in this ap- 
proximation. To correct this, we have to add satellite modes 
to the fundamental solution ( |j"8| ) 

H x (y,z) = H(y)e ikz +Hi(y,z), (19) 

z) = B+(y)e^ k+K > + B-(y)e^^ z . (20) 

The dependence of S ± versus y can be found from the 
Helmholtz equation, 



(21) 



where are constants. In order for S ± to exponentially 
decay in y, we have to assume here that k < k/2. 

Substituting Tt\ terms into the boundary condition 
( IToj ) generates first order terms that have x-dependence 
expz(/c ± k)x, and second order terms proportional to 
exp(ikx). From the former one finds that 



iknHoh 
2\/k 2 ± 2kh 



(22) 



and the latter gives an expression for the tangential electric 
field on the surface, 



1 2 Vk 2 + 2kn + Vk 2 - 2kK 

E z = —ikh kH x ■ . (23) 

4 Vk 2 - 4fc 2 

One can now solve Maxwell's equations with the bound- 
ary condition given by Eq. (23 ) (see details in [|To||). It turns 
out, that in the region of frequencies k < k/2 there exists 
a solution corresponding to a wave propagating with the 
phase frequency equal to the speed of light. The frequency 
and the loss factor of the mode are shown in Fig. [| (solid 
lines). We see that decreasing the height of the corrugation 
results in smaller wakes, and hence leads to the suppression 
of the interaction of the synchronous wave with the beam. 
In the limit of small frequencies, k <C K the frequency is 
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Figure 4: Frequency and loss factor as function of height. 

We should mention here that the perturbation theory 
breaks down for very small values of h. Indeed, we im- 
plicitly assumed that the satellite harmonics in Eq. ( |2l| ) 
are localized near the surface, otherwise our approxima- 
tion of plane surface becomes invalid. Hence, we have to 
require that k — 2k ^> a -1 , which gives the following con- 
dition of applicability: h > a -1 / 4 ^ 5 / 4 . This condition 
explains why this mode was not found by Papiernik: be- 
ing perturbative in parameter h the approach developed in 
his paper is applicable only when h can be made arbitrar- 
ily small. Let us also mention that our approach can be 
applied to the rectangular corrugations problem: we take 
Faraday's law, to obtain the averaged longitudinal field at 
r = a, (E z ) = ikSgH x /p, and then follow the method of 
Ref. to reproduce Eq. [TTJ. 

Finally, in Fig. [| we include also the results of Fig. ||, 
obtained by field matching (the dashes, the dotdashes). 
For the comparison we make the correspondences with the 
first Fourier component of the wall shape: k — 2n/p and 
h = 25 /it. We note that even though the rectangular corru- 
gations violate our requirement for smoothness, the results 
for the two methods are very similar. 
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